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ABSTRACT 

We present new observations on PKS 2155-304 with the Chandra Low Energy Transmission Grating Spec- 
trometer (LETG), using the Advanced CCD Imaging Spectrometer (ACIS). We confirm the detection of an 
absorption line plausibly identified as O VIII Lya from the warm-hot intergalactic medium associated with a 
small group of galaxies along the line of sight, as originally reported by Fang et al. 2002 (here after FANG02). 
Combining the previous observations in FANG02 and five new, long observations on the same target, we in- 
crease the total exposure time by a factor of three, and the total counts per resolution element by a factor of 
five. The measured line equivalent width is smaller than that observed in FANG02, but still consistent at 90% 
confidence. We also analyze the XMM-Newton observations on the same target, as well as observations us- 
ing the Chandra LETG and the High Resolution Camera (HRC) combination. These observations have been 
used to challenge our reported detection. While no line is seen in either the XMM-Newton and the Chandra 
LETG+HRC data, we find that our result is consistent with the upper limits from both data sets. We attribute the 
non-detection to (1) higher quality of the Chandra LETG+ACIS spectrum, and (2) the rather extended wings 
of the line spread functions of both the XMM RGS and the Chandra LETG+HRC. We discuss the implication 
of our observation on the temperature and density of the absorber. We also confirm the detection of z ~ 
O VII absorption and, comparing with previous Chandra analysis, we obtain much tighter constraints on the 
line properties. 

Subject headings: intergalactic medium — quasars: absorption lines — X-rays: galaxies — large-scale struc- 
ture of universe — methods: data analysis 



1. INTRODUCTION 

Cosmological hydrodynamic simulations predict that a 
large amount of baryons in the local universe reside in the 
intergalactic medium (IGM) (see, e.g., Cen & Ostriker 1999; 
Dave et al. 2001; Kravtsov et al. 2002; Cen & Ostriker 2006). 
While the warm and photo-ionized Lya clouds contain a 
significant fraction of the total baryons (see, e.g., Morris 
et al. 1991; Bahcall et al. 1991;Stocke et al. 1995;Shull et 
al. 1996;Penton et al. 2000; Williger et al. 2006;Lehner et 
al. 2007), the remainder is shock-heated to temperatures be- 
tween 10 5 - 10 7 K with moderate overdensities. Highly 
ionized metals in this Warm-Hot Intergalactic Medium, or 
"WHIM", can produce detectable absorption signatures in the 
UV/X-ray spectra of background sources, just as the cooler 
IGM imprints the Lya forest (see, e.g., Perna & Loeb 1998; 
Hellsten 1998; Fang, Bryan, & Canizares 2002; Chen et 
al. 2002; Viel et al. 2003; Fujimoto et al. 2004; Furlanetto 
et al.2005; Cen & Fang 2006;Kawahara et al. 2006). De- 
tectability of the WHIM gas through emission is also under 
investigation (see, e.g., Croft et al. 2001; Phillips et al. 2001; 
Kuntz & Snowden et al. 2001; Zappacosta et al. 2002, 
2005,2007; Yoshikawa et al. 2003, 2004; Kaastra 2004; Mit- 
taz et al. 2004; Slotan et al. 2005;Fang et al. 2005; Ursino & 
Galeazzi 2006;Cen & Fang 2006; Takei et al. 2007a,2007b; 
Mannucci et al. 2007). 

Recent detections of the UV absorption lines with the STIS 
onboard the Hubble Space Telescope and the Far Ultraviolet 
Spectrometer (FUSE) have firmly revealed the existence of 
the WHIM gas at temperatures between 10 5 - 10 5 5 K (see, 
e.g., Savage et al. 1998; Tripp et al. 2000; Oegerle et al. 2000; 
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Sembach et al. 2004; Tumlinson et al. 2005; Danforth & 
Shull 2005; Stocke et al. 2006). The O VI doublet at 1032 
and 1039 A can be used to probe the temperatures, densi- 
ties and redshifts of that component of the intervening WHIM 
gas. Recently, a population of intervening broad HI Lya ab- 
sorbers (BLAs) have also been detected in the far-UV band, 
stirring interests in probing WHIM with BLAs (see. e.g., 
Richter et al. 2005; 2006). The contribution of both BLAs and 
O VI absorbers to the baryonic density are ri^(BLA) > 0.0027 
(Richter et al. 2005) and fJi(OVI) ~ 0.0022 (see, e.g., Dan- 
forth & Shull 2005), respectively, assuming a Hubble constant 
of 70 km s _1 Mpc _1 . However, as simulations suggest, about 
only one third of the WHIM gas can be detected in the UV 
band — with higher temperatures the remaining two thirds 
can only be revealed in the X-ray (see, e.g., Fang, Bryan, & 
Canizares 2002; Chen et al. 2002; Viel et al. 2003; Cen & 
Fang 2006; Cen & Ostriker 2006). 

Several lines of evidence indicate that the narrow X- 
ray absorption lines from the WHIM gas may be detected 
with high resolution spectrometers onboard Chandra and 
XMM-Newton. A number of z ~ absorption lines were 
detected in the spectra of background AGNs with Chan- 
dra and XMM-Newton (Nicastro et al. 2002; Fang, Sem- 
bach, & Canizares 2003; Rasmussen et al. 2003; Kaspi et 
al. 2002;McKernan et al. 2004,2005; Cagnoni et al. 2004; 
Williams et al. 2005, 2006a, 2006b). However, current instru- 
mental resolution cannot distinguish between Galactic and the 
Local Group origin of the absorbing gas. While a number of 
studies proposed these absorption lines are produced by the 
intragroup medium in the Local Group (see, e.g., Nicastro et 
al. 2002; Williams et al. 2005, 2006a,2006b), other evidence 
suggests a Galactic-origin of this hot gas (Wang et al. 2005; 
Yao & Wang 2005; Fang et al. 2006; Yao & Wang 2007; 
Bregman & Lloyd-Davies 2007). 

Fang et al. (2002, hereafter FANG02) reported the first 
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detection of an intervening X-ray absorber at z = 0.0554 
along the sightline toward PKS 2155-304, using the Chan- 
dra LETGS. At the same redshift a small group of galaxies 
had been detected by Shull et al. (1998), who subsequently 
reported the detection of O VI (Shull et al. 2003). Subsequent 
observations with XMM-Newton did not detect the O VIII 
absorption, setting a 3a upper limit equivalent width (EW) 
of ~ 14 mA compared to our reported EW of 14. 0^ 1 mA 
(Cagnoni et al. 2004). Several O VII and O VIII absorption 
lines were detected along the sightline toward H 1821+621 
(Mathur et al. 2003), although the statistics is low and the de- 
tections are at 2 - 3a level. Mckernan et al. (2003) claimed 
detection of an intervening absorption system at z ~ 0.0147 
toward 3C 120 at > 3a level, although they cannot rule out 
the possibility that this system is intrinsic to the jet of 3C 120. 
Recently, during an extremely bright state of Mkn 421, Nicas- 
tro et al. (2005) detected two intervening absorption systems 
(z > 0) with high significance. However, observations with 
the XMM-Newton cannot confirm the detections, and the con- 
sistency between the Chandra and XMM-Newton was inves- 
tigated (Ravasio et al. 2005; Williams et al. 2006c; Kasstra et 
al. 2007; Rasmussen et al. 2007). 

In this paper, we confirm detection of the intervening ab- 
sorber toward PKS 2155-304 using new observations to re- 
examine the results from FANG02. PKS 2155-304 has been 
used repeatedly as a Chandra and XMM-Newton calibration 
target because it is one of the brightest extragalactic soft X- 
ray sources, and also because of its relatively simple spec- 
trum shape. Since we published FANG02, five more Chan- 
dra observations with the same instrumental configuration 
(LETG+ACIS) were conducted, which more than triples the 
total exposure time that was reported in FANG02. The in- 
crease in the number of detected photons (by a factor of ~ 5) 
improves the statistics of the spectrum significantly. We also 
analyze the XMM-Newton observations on the same target, 
as well as observations using the Chandra LETG+HRC com- 
bination. These observations have been used to challenge our 
reported detection. 

2. DATA ANALYSIS 

At z = 0.112, PKS 2155-304 (a = 21 /, 58 m 52.P, S = 
-30°13'32.1") is one of the brightest extragalactic X-ray 
sources. This source has been chosen to calibrate various 
instruments onboard Chandra and so is observed repeatedly. 
For our purpose, we choose observations that have been con- 
ducted with the Low Energy Transmission Gratings (LETG) 
and ACIS-S as the focal plane detector 3 , to be consistent with 
FANG02. In section §4, we will discuss results from observa- 
tions using LETG with HRC-S as the focal plane detector. 

In FANG02, we report the results from three observations 
that were conducted during May and December of 2000 and 
November of 2001, with a total exposure time of 86.7 ksec. 
Subsequently, seven more observations were conducted be- 
tween June 2002 and September 2005. We select those obser- 
vations with at least 30 ksec exposure time, which results in a 
total of five more observations. Including these five more ob- 
servations increases the total exposure time to ~ 277 ksec, and 
the number of counts in the spectral regions of interests is also 
increased by a factor of ~ 5. Table 1 lists the observation log. 
Seven of the eight observations have the nominal offset point- 
ing for LETG+ACIS-S: the aimpoint is moved by Ay = +1 .5' 
along the observatory y direction so that the most interested 

3 For LETG and ACIS-S, see http://asc.harvard.edu 



TABLE 1 
Observation Log 



Observation ID Observation Date Duration (ksec) Reference 



1703 


31 May 2000 


26.7 


1,2 


2335 


06 Dec 2000 


30.0 


1.2 


3168 


31 Nov 2001 


30.0 


1,2 


3669 


1 1 Jun 2002 


50.0 


2 


3707 


30 Nov 2002 


30.0 


2 


4416 


16 Dec 2003 


50.0 


2 


6090 


25 May 2005 


30.0 


2 


6091 


19 Sep 2005 


30.0 


2 



REFERENCES. — 1. Fang et al. (2002); 2. this paper. 



wavelength range (0 - 26 A) can be covered entirely by back- 
side chip S3, which has higher quantum efficiency and large 
effective area, and the Scientific Instrument Module (SIM) is 
moved by (SIM-Z = -8mm), to mitigate the CTI (Charge 
Transfer Inefficiency)-induced energy resolution degradation. 
This also avoids putting the zeroth-order at one of the node 
boundaries. The only exception is Obs.#3669, which had a 
y-offset of 3.3'. This put the aimpoint on chip S2. 

Data analysis is performed as described in FANG02. We 
briefly summarize here and refer the reader to that paper for 
more details. The Chandra LETGS produces a zeroth or- 
der image at the aim-point on the focal plane detector, the 
ACIS-S array, with higher order spectra dispersed to either 
side. The LETGS provides nearly constant spectral resolu- 
tion ( AA = 0.05) through the entire bandpass (0.3-5 keV). The 
moderate energy resolution of the CCD detector ACIS-S is 
used to separate the overlapping orders of the dispersed spec- 
trum. We add the plus and minus sides to obtain the first order 
spectrum. All the five data sets are analyzed with the standard 
Chandra Interactive Analysis of Observations (CIAO) 4 and 
customized Interactive Data Language (IDL) routines. 

3. CONTINUUM SPECTRAL ANALYSIS 

The continua are typically fitted with a power law modi- 
fied by neutral hydrogen absorption in the range between 6 
and 42 A. For observations #1703, 2335, and 3168 we re- 
fer to FANG02 for detailed model parameters. For #3669, 
the data is best fitted with a broken power law. The pho- 
ton indices are 2.52 ± 0.02 and 1 .96 ± 0.02, with a break en- 
ergy of ~ 2 keV, and a flux of 1.66 x 10~ 10 ergscm" 2 s _1 be- 
tween 0.5 and 2.4 keV. The remaining four observations can 
be best fitted by a single power law with photon indices of 

(2.69±0.02,2.69±0.02,2.58±0.01,2.66±0.01), for #3707, 
4416, 6090, and 6091, respectively. The flux between 0.5 and 
2.4 keV is (5.06,5.69, 10.33,6.54) x 10" 11 ergscrrrV 1 , re- 
spectively. All the models include a Galactic absorption fixed 
at Nh = 1.36 x 10 20 cm" 2 (Lockman & Savage 1995; errors 
are quoted at 90% confidence). 

In FigureQ]we show the raw counts of the five observations. 
The bin size is 0.025 A. The line spread function (LSF) of 
LETGS has a typical full width of half maximum (FWHM) of 
~ 0.05 . In this way one FWHM will contain two bins. Two 
vertical dashed lines indicate the wavelength of the redshifted 
O VIII line (at ~ 20 A) reported by FANG02, and rest O VII 

4 See http://asc.hai-vard.edu/ciao 
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FIG . 1 . — Raw counts for the five observations between 1 8 and 22 A. Bin size is 0.025 A. Two vertical dashed lines indicated positions of the redshifted O VIII 
and rest O VII lines from FANG02. 



lines (at ~ 21 .6 A). There are no strong absorption lines in 
these observations, but we do find small dips in each of these 
observations at the corresponding positions. 

To enhance the signal we stack all the eight data sets. To en- 
sure that all the eight observations are properly aligned with 
each other, we also compare the wavelengths of two known 
features, atomic O I from our Galaxy and solid-state O I from 
the instruments, and calibrate the five data sets against these 
two lines (FANG02). Only two observations (1703 and 2335) 
showed a systematic shift of ~ 0.05 A shift, and we have to 
adjust the shift manually. Figure|2]shows the stacked data be- 
tween 18 and 22 A. We have achieved from ~ 450 to 600 
counts per bin here, compared to ~ 100 counts per bin in 
FANG02. We again mark the positions of the two absorption 
lines. 

Measurement of narrow absorption features depends cru- 
cially on determining the correct continuum level. However, a 
simple power law cannot provide the best fit to the continuum 
of the stacked data for: 1) there are spectral variations among 
different observations, and 2) there are residual uncertainties 
in instrumental efficiencies. Instead we fit the local continuum 



using polynomials. More specifically, we first fit the stacked 
data with a single power law plus neutral hydrogen absorption 
as an initial guess. The residuals between 6 and 42 A are then 
fitted with a six-order polynomial. Such a polynomial will 
get rid of any features with a characteristic scale of > 6 A, but 
smaller features will be preserved. Furthermore, to account 
for small scale fluctuations, we divide the regions between 6 
and 42 A into 9 bands, with a bandwidth of 4 A each. For 
each region we then fit the residuals from the previous power 
law fitting with a polynomial. The order of this polynomial 
is given by the one that has the smallest \ 2 - The red lines in 
Figure|2]show the final result of the fitted continuum in the 18 
- 22 A range. The top panel shows the spectrum in units of 
counts per bin, with a binsize of 25 mA, and the bottom panel 
shows the flux, in units of photons cm" 2 s" 1_1 . 

How well does our model depict the observed continuum? 
We first calculate the ^-distribution of the data. \ is defined 
as: 

data count — model count 
^ (model countY I 2 

We then run 10,000 Monte-Carlo simulations, each Simula- 
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FIG .2. — The top panel shows the raw counts of the stacked data, with a binsize of 0.025 A. The bottom panel shows the flux spectrum. The red line show the 
best fitted model. Two vertical green dashed lines indicate positions of the redshifted O VIII and rest O VII lines from FANG02. 



tion is a random realization of the model based on Poisson 
statistics, and we also calculate the x-distribution of these 
simulations. Figure|3]shows the comparison between the data 
(dark line) and the simulation (red line). We find in general 
the x-distribution of the data follows that of the simulation 
reasonably. The data does show a long, negative tail, this is 
actually caused by the absorption features that we will dis- 
cuss below. Excluding the three prominent absorption fea- 
tures that we will discuss in the next section, we perform a 
Kolmogorov-Smirnov test on the two distributions. The null 
hypothesis that the data sets are drawn from the same distri- 
bution has a significant level of ~ 0.48, suggesting that our 
model fits the continuum quite well, and that our fitting pro- 
cedure does take care of small scale fluctuations. 

4. RESULTS 

A blind search of the entire region from 6 to 42 A for fea- 
tures at > 3a level give only three prominent absorption fea- 
tures, at ~ 20.0, 21.6, and 23.5 A. The last one is an interstel- 
lar O I line and will not be discussed here. These two features 
at 20.0 and 21.6 A are subsequently analyzed with customized 
routines from the software package ISIS (Interactive Spectral 
Interpretation System, see Houck & Denicola 2000) 5 . Specif- 
ically, we subtract the data from the best-fit model, and then 
fit the residual with a Gaussian line profile. The EW is an in- 
tegration of (l -///o), where / is the observed spectrum, and 
Iq is the continuum. We list the fitting parameters in Table 2. 
The Gaussian line fitting gives a line width (er = 0.017 and 

5 see |http://space.mit.edu/ASC/ISIS/ 1 



0.021 for the 20.0 and 21.6 A lines, respectively) that is even 
narrow than that of the instrumental LSF, we only present 
the 90% upper limits on the intrinsic line widths. In calcu- 
lating SNR (signal-to-noise ratio), signal is the total photon 
counts missed in the line, and noise is square-root of contin- 
uum photon counts within the line profile. The continuum 
photon counts are the total of counts under continuum, with a 
width of three times the measured line width, slightly larger 
than the FWHM of the fitted Gaussian line. In Figure [4] we 
show the residual spectra of both features (black lines) and 
the fitted Gaussian lines (red lines). 

4.1. Contamination from the Telescope Features? 

Before we proceed to further discussion, we need to under- 
stand the effects of any fixed-position detector features (such 
as those chip gaps and node boundaries) and instrumental fea- 
tures (such as those absorption edges in the telescope materi- 
als). While over the entire wavelength range the telescope has 
been well calibrated, it is still likely that some of the variations 
of the effective area and/or detector quantum efficiency across 
these features can affect the estimated significances of our de- 
tected lines, if it happens that the detected lines are close to 
these telescope features. 

We do not find any instrumental features near our two de- 
tected lines (POG 6 , Table 9.4). The positions of detector 
features depend on the aimpoint of each observation. For 
commonly used offset pointing (Ay = 1.5' and SIM-Z = 
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TABLE 2 
Line Fitting Parameters 



O viii Lya a 



O vii Hea a 



Wavelength 


20.00 ±0.01 


21.60±0.01 


cz(km s" 1 ) 


16307 ±158 


-26 ±138 


Line Width b 


< 0.027 


< 0.029 


Line Flux c 


9 1+0.6 
^• A -0.8 


3.9!?;? 


EW(m) 


7.42!?$ 


1 ~X 7<+4.09 
u - ' J -2.80 


SNR d 


5.0 


7.1 



a. Rest-frame wavelengths for O viii Lya and O vn 
Hea are 18.9689 A and 21.6019 A respectively 
(Verner et al. 1996). Errors are quoted at 90% 
confidence level hereafter. 

b. 90% upper limit of the intrinsic line width a, in 
units of A. The FWHM of the fitted Gaussian line 
is 2.35o-. 

c. Absorbed line flux in units of 
10" 5 photons cm" 2 s _1 . 

d. In units of equivalent sigma of a Gaussian 
distribution with the same confidence level as 
the Poisson significance gives, based on the x~ 
distribution. 



-8mm), the nominal aimpoint in S3 chip is moved from 
(chipx, chipy) = (23 1 , 502) to (chipx, chipy) = (49, 168), which 
also avoids the zeroth-order being at one of the node bound- 
aries. A careful calculation indicates no telescope features 
around the 21.6 A line, but such configuration does put the 
positive first order 20.0 A line close to the boundary between 
node 2 and node 3 in S3 chip. However, we believe this is 
unlikely to have large impact on the significance of the 20.0 
A line because of the following considerations. First, as we 
stated before, Obs.#3669, the one with the most photon counts 
(accounts for almost half of the total counts at 20.0 A), had a 
y-offset of 3.3', which put the aimpoint in chip S2. This re- 
sults in the positive first order 20.0 A line being far away from 
any telescope features. Secondly, this line is clearly visible in 
both positive and negative first order spectra. It is unlikely that 
this line would be detected in the negative order if it is caused 
by detector feature in the positive order. Thirdly, and most im- 
portantly, we believe that our continuum-extraction technique 
can take care of the small variation caused by node bound- 
aries. During most observations, the spacecraft is dithering 
to (1) smear out the reduced exposure from chip gaps, and (2) 
smooth out pixel-to-pixel variation in the response (see POG). 
The standard dithering pattern has a peak-to-peak width of 
16", this roughly corresponds to a response variation on the 
scale of ~ 1 for LETGS. However, based on the fitting pro- 
cedures we described in section §3, the residual between 18 
and 22 A , after broadband continuum-fitting with a power 
law and 6-order polynomial, was fitted with another 6-order 
polynomial. This would get rid of any feature with a typical 
scale of > 0.7 , which includes the dithered node boundary. 
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FIG. 3. — x-distribution of the data (dark line) and of 10,000 Monte-Carlo 
simulations (red line). 



The effect of such procedure can be viewed from Figure|2l in 
the top panel, the raw counts spectrum shows some variation 
between 19 and 21 A , however, in the bottom panel, such 
variation has been smoothed out in the flux spectrum. 

4.2. Comparison with results from FANG02 

Most importantly, our work does detect an absorption line 
at 5a significance at ~ 20 A, confirming the initial detection 
by FANG02. We find that the results presented in this paper 
are consistent with those of FANG02, but the estimated EW in 
this work (7.42!j gjj) is about one-half of that estimated earlier, 
based on many fewer counts (150 vs. 550 counts per bin) and 
with larger uncertainty. The two results are consistent at 90% 
confidence level. The properties of the stronger local {z ks 0) 
O VII absorption line are fully consistent with those reported 
in FANG02. 

4.3. Comparison with results from Chandra LETG+HRC-S 

For a consistency check, we also analyzed the observations 
up to a date of 2004 November 22 taken with the LETG oper- 
ated with High Resolution Camera for spectroscopy (HRC- 
S; Table 3). We reprocessed these observations using the 
Chandra Interactive Analysis of Observations (CIAO) soft- 
ware (ver. 3.3.0.1) with the calibration database CALDB 
(ver. 3.2.1). Because the HRC itself does not have intrin- 
sic energy resolution, the photons located to the same po- 
sition of the grating arms but from different grating orders 
cannot be sorted out as in those observations with ACIS. We 
therefore used a forward approach to account for the spectral 
order overlapping. For each observation, we first calculated 
the response matrix files (RMFs) and the auxiliary response 
functions (ARFs) for grating orders from the first to sixth 
(assuming the contributions to the spectrum from the higher 
> 6 orders are negligible). We then added these six RMF 
and ARF pairs to form an order-combined response (RSP) file 
using IDL routines from PINTofALE package 7 . To enhance 
the counting statistics, spectra from all the HRC-S observa- 
tions, after a consistency check, were co-added, and the RSP 

7 http://hea-www.harvard.edu/PINTofALE 
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FIG. 4. — Residual flux after subtracting the best-fitted continuum from the stacked data. The red lines show the ISIS-fitted Gaussian line profiles. The top 
panel is for the redshifted O VIII line and the bottom panel is for the zero-redshift O VII line. 



TABLE 3 
HRC Observations 



ObsID 


Date 




Exposure (ks) 


331 


1999 Dec. 


25 


63.2 


1013 


2001 Apr. 


06 


26.8 


1704 


2000 May 


31 


26.0 


3166 


2001 Nov. 


30 


30.0 


3709 


2002 Nov. 


30 


13.8 


4406 


2002 Nov. 


30 


13.9 


5172 


2004 Nov. 


22 


27.2 



files are weight-averaged based on the corresponding expo- 
sure time and the continuum intensity obtained from a global 
fit to each individual spectrum. 

In Fig. [5] we show the total HRC spectrum (from the first to 
sixth order) between 18 and 22 A. The 21.6 A line is clearly 
visible, but not the 20.0 A line. There are some hints of 
absorption features between 19 and 20 A, but nothing near 
the LETG+ACIS detected 20.0 A line. We cannot extract 
the first order spectrum because of the reason we described 
above, but using a first order response matrix, we estimate 
the continuum flux of the first order is about 90% of the total 
flux, which gives ~ 450 counts per 0.025 A at ~ 20.0 A re- 
gion. The observations with the LETG+ACIS have a slightly 
stronger continuum (~ 550 counts per 0.025 A at the same 



wavelength). We estimate the upper limit of an absorption 
line EW, by adding such line in the spectrum and calculat- 
ing A\ 2 by varying line EW. This gives an upper limit of ~ 9 
mA, which is consistent with the result from the LETG+ACIS 
observations. 

To further investigate the difference between the results 
from the LETG+ACIS and LETG+HRC, we compare their 
LSFs. In Figure |6] we show the LSFs of the LETG+ACIS 
(red line) and the LETG+HRC (green line) at 20 A. The blue 
line is the LSF of the XMM RGS 1 , and we will discuss it later. 
Clearly, the HRC curve has a rather broad wing. This suggests 
that when compared with that of the ACIS, more photons will 
be distributed in the wing, and the significance at the line cen- 
ter will be lower. Such broader profile is likely caused by the 
non-linearities of the LETG+HRC-S dispersion relation . To 
test the impact of the LSF, we run 10,00 Monte-Carlo simu- 
lations for both the LETG+ACIS and the LETG+HRC, based 
on the Poisson statistics. Each simulation is a random realiza- 
tion of the corresponded continuum model plus an absorption 
line at 20 A. We adopt line parameters from the O VIII line in 
Table 2 with an intrinsic line width of 0.016 A, and then fold 
it through the corresponded LSF. We find that out of 10,000 
trials, the line is detected at 3er level or higher in more than 
half of the ACIS simulated spectra (~ 57%); however, for 
the HRC spectra, the detection probability decreases to about 
one-fourth (~ 25%). 

8 For detailed discussion on the non-linearities of the LETG+HRC-S dis- 
persion relation, see http://cxc.harvard.edu/cal/Letg/Corrlam/ 
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FIG. 5 . — LETG+HRC spectrum. The red line in the top panel is the model. 
The bottom panel shows \ 2 with a sl g n - 

4.4. Comparison with results from XMM-Newton 

PKS 2155-304 was also extensively observed with XMM- 
Newton. A detailed spectroscopy study of the RGS data was 
presented by Cagnoni et al. (2004). They confirm the detec- 
tion of zero-redshift O VII absorption line at a significance 
level of 4.5ct. Their best fit equivalent width (EW = 19.50^ 
mA) is higher than what we measure but consistent at 90% 
confidence level. They also claim that this line profile is pos- 
sibly double peaked. Although the Chandra LETG-ACIS has 
higher energy resolving power, we cannot confirm this struc- 
ture in our data. 

Cagnoni et al. (2004) stated that they do not detect the red- 
shifted O VIII absorption line at ~ 20 A, reported by FANG02, 
and set a 3er upper limit of 14 mA. However, this value is large 
enough to be consistent with both the improved measurement 
and the original best fit value of FANG02 for the detected line. 

To further compare the Chandra and XMM-Newton data, 
we reanalyze the XMM-Newton data presented in Cagnoni 
et al. (2004) with the same techniques adopted in this paper. 
In doing this we hope we can minimize the differences that 
can be caused by various data analysis techniques. For in- 
stance, whereas we apply a adapted-polynomial plus power 
law fit technique to subtract continuum, they use local ab- 
sorbed power law to fit individual regions with bandwidths 
of 2-3 A. 

We use the three data sets that were analyzed in Cagnoni 
et al. (2004) (ObsID #0080940101, 0080940301, and 
0080940401). We do not select the fourth observation (ObsID 
#0080940501) because it contains a background flare. These 
three observations have a total exposure time of ~ WQksec 
(see Table 1 of Cagnoni et al.(2004) for detailed observation 
log). We analyze the data with the standard Science Analy- 
sis System (SAS) version 5.4.1 9 . We use only RGS-1 data 
since there is no data between 20 and 24 A from RGS-2 due 
to the failure of a CCD chip. The LSF of the RGS-1 has 
a FWHM of ~ 0.06 A with rather extended wings 10 . We 
will not present our data extraction procedures here, but refer 
reader to Cagnoni et al. (2004) for details. After obtaining the 
first order spectrum, we follow procedures that are described 
in section §3 to subtract the continuum. Figure [7] shows the 
RGS-1 spectrum between 19.7 and 20.7 A. We obtain ~ 450 
counts per 0.025 A bin. This is 20% lower than what we ob- 

9 See http://xmm.vilspa.esa.es/ 

10 See XMM-Newton Users' Handbook at http://xmm.vilspa.esa.es/ 
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FIG. 6. — The line spread functions of the Chandra LETG+ACIS (red), the 
Chandra LETG+HRC (green), and the XMM RGS 1 (blue) at 20 A. 

served in Chandra data (~ 550 counts per 0.025 A ). 

The dark solid line in figure [7] shows the XMM data be- 
tween 19.5 and 20.5 A. No absorption line is seen around ~ 
20 A. With our continuum subtraction techniques, we obtain 
a 3<T upper limit of ~ 10 mA, consistent with what is obtained 
in Cagnoni et al. (2004) and the ~ 7 mA equivalent width line 
we report in this paper. To demonstrate this, in Figure [7] we 
plot a 3a, 10 mA line at ~ 20 A in green. An illustration 
(but without imposing statistical fluctuations) of how the line 
detected in Chandra would appear in the XMM data is plot- 
ted in red. We adopt our Chandra line parameters, and then 
convolve the line with the RGS1 LSF 11 . 

In Figure [8] we show the A\ 2 of fitting the XMM data with 
a continuum plus a Gaussian line model (again, convolved by 
the LSF), where A% 2 = is defined as no Gaussian line. We 
fix the line center at A = 20 A, and gradually vary the line 
equivalent width. The horizontal dashed line indicates the 
Ax 2 = 9, or 3a confidence. Three vertical dark lines show 
the line width and 90% lower and upper limits of our detec- 
tion. Again, we find that one would not have expected XMM- 
Newton to detect the line seen with Chandra. 

Two main reasons that explain the non-detection of the X- 
ray absorption line in the Chandra LETG+HRC data may also 
apply there. First, the continuum is slightly higher in Chandra 
data. At 20 A, the Chandra observations show ~ 550 counts 
per bin, compared with ~ 450 counts in the XMM data. Sec- 
ondly, and most importantly, the LSF of the RGS1 not only 
is larger than that of the LETGS, but also has rather extended 
wings. As pointed out by Williams et al. (2006) in a study of 
the XMM data of Mrk 421, the central 0. 1 A region of a line 
contains about 96% of the total line flux for Chandra LETG, 
for RGS this number decreases to ~ 68%. In Figure|6]the dot- 
ted line is the LSF of the RGS1. The wings of the RGS1 are 
even broader than that of the LETG+HRC. We ran a similar 
Monte-Carlo simulation on the RGS1 data. The result sup- 
ports our finding: out of 10,000 trials, a 3cr line is detected 
in only ~ 10% RGS1 spectra, compared with nearly ^57% 
detection in the LETG+ACIS spectra, and ~ 25% detection 
in the LETG+HRC spectra. 

5. DISCUSSION 
1 1 A template of the LSF was kindly provided by A. Rasmussen. 
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FIG. 7. — XMM data at ~ 20 A. The dark line shows the data, the green 
line shows the continuum with a 3cr absorption line, and the red line shows 
the continuum with the line detected in Chandra data. The line at ~ 19.75 A 
is an instrumental feature. 

Following FANG02 and assuming that the line is not sat- 
urated, the column density of the redshifted O VIII absorp- 
tion line at 20 A is N(0 VIII)= 5.0!| % x 10 I5 cm" 2 . The line 
width sets an upper limit to the path length of ~ 5.8/170 Mpc 
based on the Hubble flow 11 . This gives «/, > 5 x 10~ 6 cm" 3 
^oi/aMjs where Z0.1 is the metallicity in units of 0.1 solar 
abundance °, /0.5 is the ionization fraction in units of 0.5 
and Z5 8 is the path length in units of 5.8/17,) Mpc. A lower 
limit to the path length of ~ 1 Mpc can be obtained by as- 
suming the absorber has the size of the small galaxy group 
detected in 21 cm images (Shull et al. 1998), which gives 
rif, w 3.1 x 10~ 5 cm" 3 Zqj/qj. This implies a range of baryon 
overdensities of 8b ~ 30- 150, which, as found in FANG02, is 
consistent with the predicted density of the WHIM gas from 
cosmological hydrodynamic simulations (see, e.g., Cen & Os- 
triker 1999; Dave et al. 2001). 

To constrain the temperature we need to study the ioniza- 
tion structure of the O VIII absorber. Assuming collisional 
ionization equilibrium, the ionization fraction is a function 
of temperature only, and /(O VIII), the ionization fraction of 
O VIII, has a peak value of ~ 0.5 between 2 and 5 x 10 6 K. 
A constraint on the temperature of the O VIII absorber can 
be obtained by studying the column density ratio between 
O VIII and O VII. Non-detection of the O VII at the corre- 
sponding redshifted position (A = 22.77 A) with a 3cr upper 
limit on the equivalent width of ~ 9.6 mA, implies N(0 VII) 
< 3 x 10 15 cm" 2 , orlog[N(0 VIII)/N(0 VII)] >0.22. This con- 
strains the temperature to be T > 2.3 x 10 6 K. However, we 
need to be careful about the assumption of collisional ioniza- 
tion equilibrium (Cen & Fang 2006). 

We can actually strengthen the constraint on temperature by 
searching for potentially co-existing absorption lines at other 
wavelengths. Recently, Shull et al. (2003) reported the de- 
tection of a series of H I and O VI absorbers along the sight 
line towards PKS 2 155-304 between 16185 and 171 16 km s _1 . 
Based on FANG02 they resorted to a multi-phase model to 

12 We use H = 70/i 7 okms _I Mpc _1 , and a standard ACDM model with 
f2 m = 0.3 and Qa = 0.7 throughout the paper 

13 We adopt solar abundance from Anders & Grevesse (1989) 



FIG. 8. — Ax 2 from fitting XMM data with a continuum plus a Gaussian 
line. The horizontal line shows the 3<r upper limit. The shadowed red area 
indicates the equivalent width and 90% confidence range of the line that is 
detected in Chandra data. 

explain the kinematic offsets among H I, O VI and O VIII ab- 
sorbers. The new measurement presented here indicates that 
these absorbers may coexist owing to the small separation in 
velocity space. For instance, the X-ray absorber may coexist 
with the component A in their observation, which is located 
between 16185 and 16252 kms" 1 and shows O VI absorption 
line. If this is indeed a single absorber, we can set a tight range 
of the OVI-to-O VIII ratio: -2.68 < log[N(0 VI)/N(0 VIII)] < 
-1.90. This in turn sets temperature: 1.6 < T < 2.3 X 10 6 K. 
Using the previous constraint from the O VHI-to-O VII ratio, 
we can put a tight constraint on the temperature of the com- 
ponent A: r«2.3x 10 6 K. 

Our discussion on the density and temperature of the 20.0 A 
absorber depends on several crucial assumptions that need to 
be closely examined. First, we assume the line is unsaturated. 
However, the 90% upper limit on the intrinsic line width can 
only give a Doppler-Z? parameter of < 566 km s -1 , much larger 
than the typical thermal broadening width ~ 50 km s~'of a 
few million-degree gas. The non-detection of the high-order 
Lyman series, particularly the Ly/3 line at the rest-frame of 
~ 16 A, can rule out the saturate scenario at la level only, and 
we notice in some cases the local X-ray absorption line could 
be saturated (Williams et al. 2006; Yao & Wang 2006, 2007). 
If the line is saturated, the derived column density can be 
treated as a lower limit only. If the line is indeed unsaturated, 
it will be broad. As suggested by the clump-infall model 
in Shull et al. (2003), a shock-wave of ~ 400 km s" 1 from 
structure evolution can provide such a non-thermal broaden- 
ing mechanism: such a shock-wave would produce a post- 
shock temperature of ~ 2 x 10 6 K. 

Secondly, we assume the gas is in collisional ionization 
equilibrium. Recent study, using numerical simulation, shows 
that non-equilibrium evolution of various ion species can 
produce significantly different effects (Cen & Fang 2006; 
Yoshikawa & Sasaki 2006), since the timescales for ionization 
and recombination are not widely separated from the Hubble 
timescale. Also, the relaxation process between electrons and 
ions may produce a two-temperature structure in the IGM and 
have impact on the metal ionization fraction, since the relax- 
ation time is comparable to the Hubble timescale (Yoshida et 
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al. 2005). However, we expect such process would be more 
important for hotter gas at T > 10 7 K (Yoshida et al. 2005). At 
very low density, especially at rit, < 10~ 5 cm" 3 , photoioniza- 
tion can also become important (see, e.g., Hellsten et al. 1998; 
Nicastro et al. 2002; Chen et al. 2003; Williams et al. 2005). 
For instance, at ri/, = 10~ 6 cm" 3 , the O VIII ionization curve 
peaks at ~ 3 x 10 5 K under the cosmic X-ray background ra- 
diation (Chen et al. 2003). 

Finally, we discuss the coexistence of O VI and O VIII ab- 
sorbers primarily because of the proximity in radial veloci- 
ties. However, as suggested by Shull et al. (2003), the O VI 
line width is narrow (40 ± 10 km s" 1 FWHM), and the sce- 
nario that O VI and O VIII can coexist is marginally plausible 
within measurement errors (logT «6.25±0.1). 

The baryonic content £!/,(0 VIII) that is probed by the O VIII 
absorption can be estimated following Lanzetta et al. (1995) 
and Tripp et al. (2000). Given the path length of the sight 
line toward PKS 2155-304 of Az « 0.116 and assuming an 
upper limit of metallicity of 0.5Z Q , we estimate 0^,(0 VIII) 
> 0.004/i7 ( ), or about 10% of the total baryon fraction. This 
number is consistent with the prediction of the WHIM gas 
from numerical simulations. Based on this single detection, 
the observed distribution dn/dz, defined as number of ab- 
sorbers per unit redshift, is higher than what simulation pre- 
dicts (see, e.g., Cen & Fang 2006, Fig. 4). We certainly need 
more detections to improve statistics. 

Note: After we submitted this paper, Williams et al. (2006b) 



published a paper on the Chandra archival data of PKS 2155- 
304. While they mainly focused on the local z = absorption 
lines, they also discussed the intervening absorption line we 
reported in FANG02. They confirmed the detection of this 
line using the Chandra LETG-ACIS. The measured equiva- 
lent width is (7.5 ±2.1 mA) is consistent with what we found, 
although the significance is lower (3.5c). They did not detect 
this line using the Chandra LETG-HRC, but obtained an up- 
per limit of ~ 12.5 mA. This is in agreement with what we 
found with the HRC data, and is consistent with the detection 
in the ACIS data. 
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